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Nomenclature 


C  specific  heat  of  the  wall .  1kg"1 

d  diameter  of  the  holes .  m 

dH  hydraulic  diameter  of  the  test  module, 

cold  side .  m 

e  wall  thickness .  m 

P[  luminosity  factor 

h  convection  coefficient .  WnrK" 

1  width  of  the  multiholed  zone .  m 

L  length  of  the  multiholed  zone .  m 


lopt  optical  path  of  the  test  module,  hot  side 

m  blowing  ratio  (ptrVtr)  l(pcVc) 

n  number  of  holes  of  the  multiholed  zone 

p  inter  row  distance 

(transverse  direction) .  m 

P  gas  pressure .  Pa 

Po  porosity  of  the  wall 

Pr  Prandtl  number 

Q  mass  flow .  kg.s"1 

r  air/fuel  mixture  mass  ratio 

s  inter  row  distance 

(longitudinal  direction) .  m 

S  heat  exchange  surface  area .  m2 

T  temperature .  °K  or  °C 

Taw  adiabatic  temperature,  hot  side .  °K  or  °C 

U  neuron  excitation 

V  velocity .  m.s"1 

UI  flux  in  infrared  units 

received  by  the  camera 
UI°  flux  in  infrared  units 

emitted  by  a  blackbody 

Greek  symbols 

OC  angle  of  inclination  of  the  holes  with 

respect  to  the  wall 

AP  pressure  loss  between  the  cold  side  and 


p  the  hot  side .  % 

At  time  step .  s 

Ax  space  between  two  pixels  of  the 

thermogram  along  the  x  axis .  m 

Ay  space  between  two  pixels  of  the 

thermogram  along  the  y  axis .  m 

e  emissivity 

<f>  radiant  flux .  W.m"2 


A,  thermal  conductivity .  W.m"1  K"1 

P  dynamic  viscosity .  kg.mV1 

p  density  of  the  wall .  kg.m’3 

°  Stefan  constant 


X  absorption  coefficient 

Subscripts 

atm  atmosphere 

c  hot  side  (pressure  =  1  atm) 

ct  casing 

ext  exterior 

f  cold  side  (temperature  =  300  K) 

hb  porthole 

jet  hole  flow 

tr  hole 

w  wall 

1.  INTRODUCTION 

The  increase  of  temperature  and  pressure  in  the  gas 
turbines  is  a  general  trend  that  allows  the  engine 
performances  to  be  increased.  These  higher 
thermodynamic  cycles  impose  cooling  systems  of 
increasingly  high  performance.  As  a  matter  of  fact,  lean 
combustion,  that  is  necessary  to  reduce  the  amount  of 
nitrogen  oxides  (Nox)  produced,  consumes  more  and 
more  air,  and  this  air  is  taken  mainly  from  the  air  usually 
reserved  for  cooling.  The  decrease  of  the  quantity  of  air 
available  for  the  cooling  systems  and  the  increase  of  the 
gas  temperature  in  the  combustion  chamber  are 
contradictory  elements  of  the  problem  which  oblige  us  to 
better  control  wall  temperatures  with  a  minimum  air 
consumption. 

Multiholed  walls  in  combustion  chambers  remain  a 
satisfactory,  efficient  and  low-cost  solution  to  the 
problem.  We  have  designed  a  test  cell  which  allows  us  to 
improve  our  knowledge  of  this  type  of  cooling  system  by 
varying  the  aerothermodynamic  and  geometric  parameters 
[1].  Many  studies  have  been  conducted  on  a  single  hole 
[2,3]  or  on  one,  two  or  three  rows  of  holes  [4,5].  These 
fairly  complete  studies  well  show  the  aerodynamic  (jet 
penetration,  turbulence  analysis)  and  thermal  (localized 
increase  of  heat  transfer  in  and  around  the  hole)  aspects, 
but  are  in  most  cases  applicable  to  the  cooling  of  turbine 
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vanes  where  the  blowing  ratio  is  low  (m  =  2)  compared 
with  those  found  in  multiholed  combustion  chamber  walls 
(m  =  3  to  5). 

Our  experimental  study  concerns  plate  shapes 
encountered  in  combustion  chambers  at  scale  1  in  a 
realistic  aerodynamic  and  thermal  environment.  Some 
studies  on  this  type  of  surface  exist  but  were  conducted  at 
different  scales  [6],  at  lower  temperature  levels  [7]  or  also 
using  a  different  mode  of  supply  of  the  holes  [8].  In  our 
case,  the  cooling  air  does  not  completely  flow  across  the 
holes  (see  figure  7),  which  corresponds  exactly  to  the  type 
of  air  flow  encountered  in  combustion  chambers.  We 
have  carried  out  a  number  of  tests  in  order  to  determine 
the  influence  of  the  parameters,  shown  in  figure  2,  on  the 
cooling  of  the  multiholed  wall.  The  use  of  neural 
networks  emerged  naturally,  compared  with  traditional 
correlation  methods,  due  to  the  multiplicity  of  tests  and 
the  large  number  of  parameters  considered.  This 
constitutes  a  specific  use  of  the  neural  networks  which  are 
usually  reserved  for  problems  such  as  command  control 
or  character  recognition. 


The  geometric  and  aerothermodynamic  parameters  under 
study  vary  as  follows: 


geometric  parameters: 

aerothermodynamic 

parameters: 

20<a<  120°; 

0.3  <  d  <  0.7  mm  ; 

1.025  <  p  <  1.685  mm  ; 

0.905  <  s  <  4.03  mm. 

0.5  <  AP/P  <  12  %  ; 

600  <  Tc  <1400  K  ; 

10  <  Qf  <  400  g/s  ; 

35  <  Qc  <  500  g/s. 

The  experimental  study  has  been  conducted  on  18  plates 
of  different  shape,  selected  so  as  to  enable  each  of  the 
geometric  parameters  (angle  and  diameter  of  the  holes, 
transverse  and  longitudinal  pitch)  to  be  studied  separately. 
For  each  plate,  we  have  studied  33  specific 
aerothermodynamic  conditions  selected  in  order  to  scan 
the  aerothermodynamic  characteristics  encountered  in  the 
combustion  chambers  of  gas  turbines. 


Primary  air  (hot) 


P 

i 

**  s 
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Figure  2.  Geometric  parameters. 


The  measurements  were  made  using  an  infrared  camera. 
Our  problem  can  be  represented  schematically,  as  shown 
in  figure  3,  by  assuming  that  we  aim  at  an  object  under  a 
temperature  T0,  in  a  given  environment,  and  this  through 
an  IR  porthole  and  a  gaseous  atmosphere. 


Figure  3.  Schematic  representation  of  the  IR 
measurement  problem. 


If  we  call  the  radiometric  flux  measured  by  the  infrared 
camera  <f>,„,  we  may  write: 


^  aim  '''  0  &(T0)  + 

^  atm  Po  ®{Te  ) 
wall  environment 

+  Tatm  TO0(Tf)+^~  T  ami  )  HTatm  ) 

v _ >■ - v - ' 

bottom  atmosphere 

+  (l-r  h)&{Th) 


window 
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On  our  thermograms,  processed  and  stored  in  the  form 
of  computer  files,  we  may  see  the  cooling  of  the  plate 
obtained  as  a  result  of  the  presence  of  the  multiple  holes. 

Since  our  aim  is  to  define  correlations  that  can  be  used  to 
presize  the  cooling  system  of  the  combustion  chamber 
walls,  we  have  tried  to  synthesize  all  the  measurement 
results  by  making  a  global  modelling  of  the  multihole 
cooling  system. 


2.  THERMAL  MODELING: 


2.1.  Flux  balance: 


The  flux  balance  of  the  multiholed  plate  applied  to  a 
pixel  of  the  infrared  image  having  the  co-ordinates  (x,  y) 
is  written  using  the  following  equation  of  energy: 

'  d2Tw(x,y,z,t) 
dx2 

d2Tw(x,y,z,t) 

T  a? 

d2Tw(x,y,z,t) 

y  dz2  J 

The  plate  thickness  is  0.8  mm.  It  may  therefore  be 
considered  as  thin:  the  temperature  gradient  along  z  is 
thus  considered  to  be  zero. 

Which  allows  us  to  write: 


dTJx,y,z,t)  _ 
v  w  dt 


P  wCwe 


dTw(x,y,0,t) 

dt 


=  4> 


cd,x 


Kd,y  +  t  cv,c 


+  <l Kay,c  ~§cv,f  -§ray,f  ~§cv,tr 

where: 

•  <sfcd  x :  density  of  heat  flow  rate  transferred  by 

conduction  along  x  through  the  wall; 

•  tycd  :  density  of  heat  flow  rate  transferred  by 

conduction  along  y  through  the  wall; 

•  (j>cv  c  :  density  of  heat  flow  rate  transferred  by 

convection  of  the  burnt  gases  towards  the  wall; 

•  <\)  ray  c :  density  of  heat  flow  rate  transferred  by 

radiation  of  the  burnt  gases  and  the  hot  casings 
towards  the  wall; 

•  c \>cvj :  density  of  heat  flow  rate  emitted  by 

convection  from  the  wall  towards  the  cooling 
gases; 

•  §rayj :  density  of  heat  flow  rate  emitted  by 

radiation  from  the  wall  towards  the  cold  casings 

and  the  porthole; 

•  <sfcvtr :  density  of  heat  flow  rate  emitted  by 

convection  towards  the  cooling  gases  flowing 
across  the  holes  of  the  multiholed  zone. 


the  holes  follows  the  Colburn  formulation  taking  into 
account  the  relative  surface  area  of  the  hole  with  respect 
to  the  plate. 


\  tr  K  U  y )  Vw  U  >’A  1 )  -  Tf  U  >’)] 

where: 


nde  V  .  .  nd 
Ax  Ay  - 


v  since. 


A-1 


4sina 


•  htr(x,  y)  according  to  Colburn; 


•  =  0  :  the  pixel  has  no  multiholed  plate  hole; 

•  L,tr  =  1 :  the  pixel  has  a  multiholed  plate  hole. 


Using  the  equations  of  the  flux  balance  written  in 
transient  mode  (the  temperature  of  the  hot  gases  rises 
briefly,  which  causes  the  plate  temperature  to  rise)  and  in 
steady-state  mode  allows  us  to  calculate  hc  and  Taw ,  the 
only  unknowns  of  the  problem,  related  by  the  relation 
<l>cv,c  =hAx,y)[Taw(x,y)-Tw(x,y,0,t)\,  according  to  a 

method  based  on  J.  Descoins’  method  [9], 

This  method  is  applied  to  all  of  our  experimental  cases, 
which  allows  us  to  know  hc  and  Taw  profiles  for  each 
aerothermodynamic  condition  corresponding  to  a  test. 


2.2.  Modeling  he  and  Taw: 

Knowing  hc  and  T;lw  ,  it  remains  for  us  to  model  them  in 
order  to  predict  wall  temperatures  using  the  flux  balance. 
We  build  this  model  from  all  our  experimental  results. 
The  relative  influence  of  hc  on  the  calculated  wall 
temperature  value  is  far  less  than  the  Taw  value.  For  this 
reason,  we  have  modeled  hc  using  a  simple  method  (see 
figure  4)  and  have  given  our  attention  to  Taw. 


he  profile  and  average 


<tw-  <bcd,y>  Pray.c  ,  <1 >cv,/  and  Pray,/  are  modeled  using  a 
conventional  method:  finite  differences,  radiation 
according  to  the  usual  formulation  and  convection 
coefficient  according  to  Colburn.  The  convection  flux  in 


Figure  4.  Experimental  profile  of  hL  and  value  modeled  in 
the  multiholed  zone. 
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The  figure  5  shows  4  longitudinal  profiles  of  convection 
coefficients,  on  the  hot  side,  obtained  from  the  IR 
measurements  and  calculated  from  the  flux  balance. 
These  profiles  were  obtained  during  4  tests  corresponding 
to  4  different  aerothermodynamic  conditions  conducted 
on  a  same  plate  and  show  how  difficult  it  is  to  model  hL. 


Figure  5.  Longitudinal  profiles  of  hc  . 


♦  For  x  >  L  (downstream  of  the  multiholed 
zone) 


fa,/ 


-n  e) 


1  -  exp 


-  L 


i  \c< 


2000  Cs 


+  11, 


Where: 

•  T\u :  adiabatic  efficiency  upstream  of  the  multiholed 
zone  ( upstream  zone).  Although  its  value  is  zero  in  theory 
(in  this  case  the  wall  is  solid  and  the  adiabatic 
temperature  is  equal  to  the  temperature  of  the  hot  gases: 
Taw  =  Tc),  we  admit  a  value  slightly  different  from  zero  in 
order  to  compensate  for  the  errors  due  to  the  convective  - 
radiative  transfer  models  used. 

•  rp,  :  adiabatic  efficiency  at  the  end  of  the  multiholed 
zone  (end  zone),  given  by  the  following  relation: 
rp,  =T\aw{x  =  L  ).  It  represents  the  maximum  adiabatic 
efficiency  value. 


We  observe  a  large  variety  of  profiles  of  convection 
coefficients  on  the  hot  side  and  we  have  therefore  chosen 
to  model  the  difference  Ah  between  the  average  value  of 
hc  calculated  in  the  multiholed  zone  (see  figure  4)  and  the 
value  of  h  calculated  according  to  Colburn  corresponding 
to  the  value  in  the  zone  located  upstream  of  the 
multiholed  zone,  as  follows: 


Ah  =  165(1 -c  ) 


where  : 


x  =  Pjet  xV.je,  xPoxL-  ; 


f  „  \ 


*  P  jet  P  / 


Pf 

V  1  J 


V .  = 

jet 


2  A P 


jet 


„  nxnxd  .  ,T 

•  Po  = -  /  (Lxl). 


Each  calculated  adiabatic  temperature  is  modeled  from 


the  adiabatic  efficiency  1]  aw 


according  to  a 


mathematical  formulation  in  which  8  coefficients  are 
used.  The  abscissa  x  and  the  length  of  the  multiholed  zone 
are  dimensionless  in  x*  and  L*  in  accordance  with 
reference  values. 

Mathematical  formulation  of  the  adiabatic  efficiency 
modeling: 

♦  For  x  <  L  (in  the  multiholed  zone) 


.(•'•')  (I  n.,)|l;-exp  - 


K  +44,604 
2000C, 


1  — 


+  44,604 
2000 


+n„ 


•  r|rf  :  adiabatic  efficiency  downstream  of  the 
multiholed  zone  ( downstream  zone)  given  by  the 
following  relation:  r|(/  =  (l  -  |Arp;  |)rp. 


For  each  test  case,  the  coefficients  Ci=1 6  as  well  as  the 
r|((  and  Arp,  values  are  calculated  using  the  Nelder 
Mead’s  method  of  the  simplex  which  calculates  the 
8  coefficients  of  the  adiabatic  efficiency  curve  whose 
corresponding  wall  temperature  curve  is  the  closest  to  that 
obtained  with  the  infrared  camera. 

The  figure  6  shows  the  relevance  of  adiabatic 
temperature  mathematical  modelling  using  the  8 
coefficients. 


X 


Figure  6.  Actual  and  modeled  adiabatic  efficiencies. 
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These  8  coefficients  depend  on  the  geometric  and 
aerothermodynamic  parameters  and  are  calculated  for 
each  test  using  Nelder  Mead’s  solver,  and  working 
therefore  in  a  space  of  8  dimensions.  The  progression  of 
the  solver  is  achieved  by  reducing  progressively  the 
difference  between  the  actual  temperature  of  the  plate  and 
the  temperature  simulated  from  the  flux  balance  where  the 
only  unknown  is  the  plate  temperature.  At  each  iteration, 
a  simulation  is  thus  calculated  from  the  corresponding 
8  coefficients. 

In  order  to  simulate  a  plate  temperature,  we  have 
therefore  to  build  a  model  with  the  8  coefficients.  Due  to 
the  complexity  of  a  mathematical  formulation  of  the 
8  coefficients  from  geometric  and  aerothermodynamic 
parameters,  we  have  moved  towards  the  use  of  neural 
networks.  In  fact,  neural  networks  are  capable  of 
interpolating  complex  functions,  without  mathematical 
formulation.  After  a  so-called  learning  phase,  we  may 
integrate  them  into  modules  of  existing  software 
programs  (in  this  case  the  program  we  have  written  for 
calculating  the  wall  temperatures),  which  allows  the 
problems  to  be  handled  in  cases  where  conventional  data 
processing  systems  are  too  weak  [10].  We  therefore  use  8 
neural  networks,  one  for  each  coefficient  to  be  simulated. 

3.  INTRODUCTION  TO  NEURAL 
NETWORKS: 

3.1.  The  formal  neuron: 


In  most  cases,  the  connections  serve  to  calculate  the  input 
of  the  neuron  whereas  the  activation  functions  determine 
the  output. 

3.3.  Me  Culloch  and  Pitts’  neuron  -  the 
multilayer  perceptron: 

The  first  definition  of  the  neuron,  that  of  Me  Culloch  and 
Pitts,  dates  from  1943. 


Therefore,  the  input  function  can  be  calculated  very 

n 

easily:  Ui  =  •  Xj  where: 

j= i 

•  Wjj  is  the  weight  of  the  interconnection  between 

the  n  upstream  neurons  j  and  the  downstream 
neuron  i; 

•  Xj  is  the  excitation  of  upstream  neuron  j; 

•  I/,  is  the  input  of  the  neuron  being  considered. 


A  neuron  can  be  described  as  a  small  decision-making 
automaton.  It  consists  of  2  parts  having  the  following 
distinct  functions: 

•  evaluation  of  the  stimuli  received; 

•  evaluation  of  its  internal  state. 


External 

signals 


Evaluation  of 

U 

Evaluation  of 

received  signal 

internal  state 

fonction  E  fonction  f 


Figure  7.  Formal  neuron  diagram. 


A  neuron  is  therefore  characterized  by: 

•  its  state  (X) ; 

•  the  excitation  level  received  at  its  input  (U) ; 

•  its  transition  function  f ; 

•  its  input  function  E  (a  sum  in  general). 


3.2.  The  interconnections: 

The  interconnection  topology  of  the  neural  networks  can 
be  extremely  varied.  Several  types  of  networks  can  be 
created  simply  by  modifying  the  interconnection  rules. 
The  singularity  of  the  interconnection  is  to  be  weighted. 
This  allows  the  relationship  between  the  neurons,  i.e.  the 
influence  of  the  activation  of  a  neuron  on  another  neuron, 
to  be  simulated. 


The  activation  of  the  neuron  is  therefore  given  simply 
by  v,  =  /({/,  ).  In  our  case,  f  is  not  the  thresholding 
function  but  the  sigmoid  function 

Xj  =  f  (U  / )  = - - - ,  which  is  typically  used  in 

l  +  e~  ' 

scenarii  like  ours  and  which  is  close  to  the  thresholding 
function. 


The  network  presented  herein,  the  multilayer  perceptron 
(see  figure  9)  which  is  the  most  appropriate  to  the  study 
of  our  problem,  has  the  following  characteristics: 

•  a  single  input  layer; 

•  a  single  output  layer; 

•  it  may  contain  one  or  more  hidden  layers; 

•  each  neuron  is  only  connected  to  all  the  neurons 

of  the  next  layer. 


Input  layer  Hidden  layers  Output  layer 

Figure  9.  Example  of  multilayer  perceptron. 
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Although  this  type  of  model  reproduces  faithfully  the 
biological  neuron,  its  implementation  implies  that  the 
weights  associated  with  the  connections  are  known.  A  so- 
called  learning  phase  must  therefore  be  applied  to  the 
network  to  set  the  weights  of  the  connections. 

3.4.  The  learning  phase: 

In  order  that  the  network  performs  the  best  interpolation, 
inputs  are  transmitted  to  the  network  and  the  network  is 
requested  to  modify  its  weights  to  obtain  the 
corresponding  outputs.  The  outputs  in  question  are 
obtained  from  the  Nelder  Mead’s  solver  applied  to  our 
test  conditions.  The  algorithm  consists  in  first  propagating 
the  inputs  forwards  until  an  output  calculated  by  the 
network  is  obtained.  The  second  step  compares  the 
calculated  output  with  the  known  actual  output.  The 
weights  are  then  modified  so  that,  when  the  next  iteration 
is  performed,  the  error  made  between  the  calculated 
output  and  the  known  output  is  minimized.  Nevertheless, 
hidden  layers  shall  not  be  forgotten.  For  this  reason,  the 
error  made  is  propagated  backwards  up  to  the  input  layer 
while  the  weights  are  being  modified.  This  process  is 
repeated  on  all  the  examples  until  an  output  error 
considered  to  be  negligible  is  obtained. 

We  have  used  the  software  program  SNNS  from  the 
University  of  Stuttgart  to  manage  all  neural  network 
issues.  The  software,  quite  user  friendly,  enables  the  user 
to  create  a  network  geometry,  perform  the  learning  phase 
(several  algorithms  are  available)  then  use  the  network  for 
function  interpolation  [11]. 

4.  USING  THE  NEURAL  NETWORKS: 

4.1.  Inputs: 

As  explained  above,  each  of  our  8  neural  networks  have 
4  inputs  corresponding  to  4  aerothermodynamic 

T 

parameters,  AP/P,  pcVc,  pfVf  and  — ,  and  one  output, 

Tf 

the  modeled  coefficient.  Each  input  and  output  is 
maintained  within  the  interval  [0,1]  in  which  the  network 
works  for  numerical  reasons  by  using  the  corresponding 
minimum  and  maximum  values  as  follows: 

_  x~  ''min 
*[0,1]  “  ~  ~  • 

-'max  -'min 

Thus,  this  enables  a  plate  of  given  shape  to  be  modeled. 
The  possibility  of  modeling  all  the  shapes  and 
aerothermodynamic  conditions  may  be  envisaged  by 


adding  geometric  parameters  to  the  inputs  of  the  neural 
networks. 

4.2.  Geometry  of  hidden  layers: 

One  of  the  crucial  problems  of  using  the  neural  networks 
in  good  conditions,  is  the  number  of  hidden  layers  as  well 
as  the  number  of  neurons  on  each  layer.  If  the  number  of 
neurons  is  too  low,  the  learning  phase  will  not  converge, 
the  network  will  have  two  many  data  to  learn  and  not 
enough  neurons  to  store  them:  the  network  will  tend  to 
underfit  the  interpolations;  it  is  not  enough  intelligent.  On 
the  contrary,  if  the  number  of  neurons  is  too  high,  the 
network  will  converge  rapidly  but  the  interpolated  points, 
excluding  the  points  learnt,  will  be  completely  false:  then 
there  will  be  overfitting:  the  network  is  too  intelligent. 

In  order  to  compare  the  effectiveness  of  the  different 
networks,  we  have  performed  so-called  propagation 
phases  which  consist,  after  the  learning  phase  has  been 
completed  on  a  given  known  function,  in  propagating 
forwards  the  input  values  across  the  network  so  as  to 
calculate  the  output  values  and  to  compare  them  with  the 
given  known  function. 


5.  RESULTS 


The  figure  10  gives  the  results  of  the  simulation  per 
neural  network  of  a  plate  temperature  profile  as  well  as 
the  corresponding  adiabatic  efficiency  on  a  case  which 
was  not  used  for  learning. 

The  plate  temperatures  measured  by  infrared 
thermography  and  the  simulation  results  are  very  close,  in 
particular  in  the  case  of  the  two  important  parameters 
which  characterize  the  cooling  system:  the  minimum 
temperature  reached  and  the  gradient  along  the  plate. 
Therefore,  we  observe  the  ability  of  our  model  to  predict 
the  temperature  of  multiholed  flat  plates  placed  between 
two  given  flows  and  whose  geometric  and  aerothermal 
parameters  are  part  of  our  field  of  study. 

Similarly,  we  observe  the  efficiency  of  the  solution  given 
by  the  neural  networks  to  a  problem  like  ours,  when  a 
mathematical  formulation  seems  difficult  to  write  in  the 
form  of  a  correlation. 

The  unknowns  of  our  problem,  hc  and  Taw  are  then 
represented  by  the  two  parameters  used  in  the  formulation 
of  Ah  and  by  the  8  neural  networks  which  provide  us  with 
an  adiabatic  efficiency  curve. 
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Figure  10.  Simulation  of  the  temperature  of  a  plate  not  used  in  the  learning  process 


On  figure  11  below,  we  present  a  typical  measurement 
and  a  comparison  with  the  calculation.  We  can  observe 
the  cooling  effect  within  the  multiholed  zone  and  its 
prolongation  after  this  zone.  The  film  cooling  created 
after  the  multiholed  zone  is  very  interesting  because 
there  is  no  air  consumption  to  protect  the  wall. 


X(mm) 


Figure  11:  Temperature  measurement  along  a 
multiholed  plate 

The  figure  12  shows  different  temperature  profile  along  a 
multiholed  plate  with  different  diameter  of  holes.  Of 
course  the  gradients  are  greater  when  the  diameter 
increase  but  the  air  consumption  is  increasing  too.  We 
can  observe  that  0.5  mm  diameter  is  a  good  compromise 
and  avoid  an  augmentation  of  the  temperature  given  by 
lower  diameter  in  the  first  rows. 


Figure  12:  diameter  influence  on  the  temperature 
profile 


The  influence  of  the  longitudinal  pitch  is  given  figure  13 
below.  When  the  pitch  decreases,  the  air  consumption 
and  the  gradient  of  the  temperature  profile  increase.  A 
combination  of  the  longitudinal  pitch  and  the  hole 
diameter  is  interesting  to  do,  in  order  to  control  the 
cooling  effect  and  the  air  consumption. 

DP-Edtl. 


Figure  13:  longitudital  pitch  influence 


OP  -2  «  ■*, 


On  figure  14,  is  given  the  influence  of  the  hole  angle.  It 
is  surprising  to  see  that  the  hole  angle  has  no  determinant 
influence  on  the  cooling  effect  within  the  multiholed 
zone.  This  means  that  the  heat  transfer  by  convection  in 
the  holes  is  not  significant  regarding  the  wall  protection. 
The  more  important  characteristic  to  protect  the 
combustion  chamber  wall  is  the  development  of  a  film 
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cooling  along  the  plate  that  preserves  this  wall  from  the 
hot  gases 


6.  CONCLUSIONS 

We  have  used  the  flux  balance  in  a  reverse  way: 
knowing  all  the  fluxes,  with  hc  modeled  by  the  relation 

Ah  =  165(1  —  e_100v  )  and  Taw  obtained  via  neural 

networks,  the  only  unknown  of  the  problem  is  the  plate 
temperature.  It  can  therefore  be  calculated. 

Thus,  a  plate  temperature  profile  can  be  simulated  for  an 
aerothermodynamic  case,  provided  that  they  are  part  of 
the  field  of  study  of  the  parameters  of  the  problem,  the 
field  of  study  in  which  the  neural  network  is  capable  of 
making  an  interpolation  of  an  adiabatic  efficiency. 

This  study  was  extended  to  the  influence  of  he  geometric 
characteristics  of  the  multihole  cooling  device  which  are 
the  hole  diameter,  the  longitudinal  pitch  and  the  hole 
angles.  We  have  measured  these  influences,  and  this 
gives  us  the  ability  to  control  the  temperature  gradient 
and  the  air  consumption  for  a  numerous  different  design 
of  multihole  cooling. 
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Paper  Number:  10 

Name  of  Discusser:  B.  Simon,  MTU  Aero  Engines  Munich 

Question: 

You  detected  no  influence  of  inclination  on  wall  temperature.  I  would  conclude  that  effusion  holes  were  to 
narrow  spaced. 

Answer: 

The  experiments  were  done  with  high  blowing  ratio  (2  to  5  or  more).  I  suppose  the  result  would  not  be  the 
same  for  lower  blowing  ratios,  but  we  observed  that  this  is  no  significant  effect  due  to  the  augmentation  of 
the  hole  surface  because  of  the  inclination.  Perhaps  there  is  a  saturation  of  the  cooling  effect  in  the  holes 
due  to  the  high  blowing  ratio. 


Name  of  Discusser:  H.  Weyer,  DLR  Cologne 

Question: 

Did  you  care  of  the  effect  of  flow  properties  on  the  “cold  side”  on  the  heat  transfer? 
Answer: 

Our  test  facility  don’t  allow  us  to  change  the  properties  of  the  cold  air 
This  cold  air  have  ambient  properties  atmospheric  pressure. 
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